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The present paper proposes a stochastic 
model of the traffic flow. This model has 
a discrete set of states and the continuous 
time. The model is a generalization of the 
discrete stochastis model that has been con- 
sidered in a previous paper of this authors 
collective, for the case of a mixed traffic flow 
which consists of "fast" and "slow" vehicles. 
Moreover the continuous time model allows 
to developed the manner for estimation of 
characteristics of the traffic flow. 

1. Introduction 

The analysis of the traffic flow and problems of the traffic con- 
trol becomes more pressing last years. Deterministic (hydrody- 
namic) models, that do not take into account some substantial 
factors as availability of several lanes, the different character- 
istics of traffic participants and so on, parallel with stochastic 
models are considered. In stochastic models the road is treated 
as a set of cells that proportional to the dynamic distance and 
a discrete time unit. The traffic flow velocity consists of the 
deterministic component connected with dynamic distance and 
a stochastic component that individual behavior into account. 
So it has been studied the particles of different types movement 
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on a set of cells. The dimension of this set is N X m, N » 1, 
m > 1. 

The discrete model of traffic flow on a one-lane road has been 
considered in [1,2]. The road is represented by a set of subse- 
quent cells. In each cell there is one or no particle (vehicle). 
The transitions of particles are possible only at discrete time 
units nA. Suppose that a cell is occupied by a particle and the 
cell following ahead is vacant. Then with probability depending 
on the particle type (vehicle velocity) transition of the particle 
forwards to the following cell occurs. The steady characteristics 
of the considered system have been investigated. The multilane 
simulation models have been considered in [3]. The analysis of 
traffic flow characteristics on the basis of simulation has been 
presented in [4, 5]. Traffic flow mathematical models have been 
investigated in [6, 7]. A deterministic model has been investi- 
gated in [6] . It has been considered a stochastic model in terms 
of the queuing theory in [7] . 

A two-lane traffic flow discrete stochastic model, that is a 
generalisation of the models, considered in [1,2], has been devel- 
oped in [8]. In this model a road is represented by two closed 
sequences of cells (lanes). The transitions of particles are pos- 
sible only at discrete time units nA. If the cell following ahead 
is busy and the two adjancent cells are empty then with the ap- 
propriate probability the particle is rebuilt in adjacent lane with 
propulsion. In [8] the approximate method has been elaborated 
for evaluation of the steady state probabilities of considered sys- 
tem fragments and such characteristics as the particle flow rate 
(traffic flow rate) and the changes of lane rate. The exactness of 
proposed approximations has been estimated by means of sim- 
ulation modelling. 

In the present paper the model is investigated which differs 
from the model considered in [8] by that in this model transitions 
of particles occur at arbitrary instants with an assigned intensity. 
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There are two types of particles characterized by the intensity 
of particle transitions. Thus the model takes into account the 
availability in the traffic flow of "fast" and "slow" vehicles. 

The continuous time model is the limit case of the discrete 
time model when A — > and accordingly the transition proba- 
bilities p — > 0. 

The method used for calculation of the studied character- 
istics is similar to the method proposed in [8]. It is proved 
that the continuous time model allows to deal with less compli- 
cated calculations than the appropriate discrete time model. It 
is explained by that in a continuous time model every change of 
states of a cells set is realized by a transition of the only particle. 

2. Model description 

Suppose that there are N x 2 cells. Two indexes correspond to 
each cell. The value of the first index can be equal to 1, 2, . . . , N 
and the second index (lane index) can be equal 1 or 2. So the 
pair of numbers i = 1, 2, . . . , N, j = 1, 2, corresponds to 

the cell. There are M\ particles of the first type and M2 particles 
of the second type. At each instant each particle occupies one 
cell and no cell contains more than one particle. 

Suppose that the value of the first index, which equal to N+l, 
is identified with 1 and the second index, which is equal to j + 1, 
is identified with 1 when j ' = 2. 

The particle movement is submitted to the next rules. If at 
time t a particle of the kth type (k = 1, 2) occupies the cell 
and the cell (i + 1, j) is vacant then during the time interval 
(t,t + A) with probability /L^ A + o( A) , A — » 0, the particle passes 
from the cell to the cell (i + 1, j) (i.e. the particle move 

forward remaining in the same lane). If the particle of the kth 
type occupies the cell the cell j) is busy, and the cells 

(i,j + 1), (i + 1, j + 1) are empty then during the time interval 
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(t, t + A) with probability //fcA + o(A) £/ie particle moves to the 
cell (i + l,j + l) ; i.e. moves forward with change of lane. In the 
other cases the particle is not replaced. The probability of two 
or more transitions during the time interval is the infinitesimal 
of higher oder than A. 

3. Evaluation of the steady state probabilities 

Let us consider the four cells with the first index, which is equal 
to 1 and 2. For approximate calculation of the steady states of 
the four cells we use method similar to the method which has 
been elaborated in [8]. 

The mean idea of the approximated method is that the consid- 
ered four cells are described by a Markov process with discrete 
set of state and continuous time. The transition intensity for 
each pair of the process is calculated provided that in every cell, 
adjancent to the four cells, ((iV, 1), (N, 2), (3,1), (3,2)) with 
probability p\ = m\/{2N), that is independent on the other cells 
states, there is a particle of the first type, with the probability 
P2 = 777-2/ (2 N) the cell is occupied by a particle of the second 
type and with probability 1 — p (p = p\ + p^) the cell is empty. 

Taking into account that the configurations of particles, sym- 
metrical respecting to line between the lanes, are characterized 
by just the same probabilities, we compound them to one state. 
There are 45 different states of the four cells shown in fig. 1. 
Suppose that equality 0(7, j) = k, k = 1,2, means that the cell 
(7, j) contains the particle of the kth type. The equality 6(i,j) = 
means that the particle (7, j) is empty. Then, for example, the 
states Ei and E2 are described as Ei = {0(1,1) = 0(1,2) = 
0(2,1) = 0(2,2) = 0}, E 2 = {0(1,1) = 0(1,2) = 0(2,1) = 
0, 0(2, 2) = 1} U{0(1, 1) = 0(1, 2) = 0(2, 2) = 0, 0(2, 1) = 1}. 

Let Pi(t) be the probability that at time the state of the four 
cells is Ei. The considered Markov process satisfies the condi- 
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Figure 1: States of the four cells 



tions of the ergodicity [9] . So the steady state probabilities pi 
exist, pi = lim PAt). 

Our assumptions allow to obtain the equations for the steady- 
state probabilities by the standard manner [9] . For example the 
equations corresponding to the states E\ and are 

-2(pi/ii + p2l^2)pi + (1 - P 2 )H\V2 + (1 - p 2 )^2Pi = 0, (1) 

-((l-p 2 )/Jl + 2(pi/J 1 +p 2 IJ2))P2+Pm+'2('L-p)PlP6+('L-p)P2P7 = 0. 

(2) 

The system for steady-states probabilities consists of the 45 
equations, corresponding to 45 system states, and the normal- 



izing equation 
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(3) 



4. Macro-characteristics of flow 

According to the Markov processes theory each equations for the 
appropriate state is a consequence of the other ones. The system 
of the 45 equations and the normalizing equation has the only 
solution i = 1, ... ,4. This solution yields the approximate 
values of the steady state probabilities of the four cells. 

As it has been mentioned above the density pk of the flow 
of particles of the kth type flow satisfies the relations pk = f^. 
k = 1, 2. On the other hand 
1 . 

Pi = 2(P4+2p9+Pl0+Pl2+Pl3+Pl6+Pl7+P20+P21+P22+P23+2p28+ 
2]?29 +^30+^31 +^32 + P33 + 2^36 + 2l?37 + 2^38+^39 +^40 + P41 +2?42) , 

(4) 

1 

P2 = 2(^5+^10+2PH+Pl4+Pl5+Pl8+Pl9+P24+P25+P26+P27+P30+ 
P31 +.P32 + ^33 + 2l?34 + 2^35 + P39 +2?40 +P41 +2?42 + 2^43 + 2^44 + 2^45 ) • 

(5) 

Let us prove the relations (4) and (5). The pair of cells (i, 1), 
(i,2) is identical with the pair (j, 1), (j, 2), 1 < i,j < AT. So 
the mean number of particles of kth type (/c = 1,2) in each 
pair equals M&/A", i.e. equals 2pk- On the other hand the mean 
number of the particles of the kth type in the pair of cells is 
equal to the sum of probabilities of the states of the considerd 
four cells in that there is exactly one particle in the pair of cells 
and the double probabilities of the states with two particles in 
this pair of cells. Hence the equations (4) and (5) are true. 

Let p k be the value of p k which satisfies the relations (4) and 
(5) if in this relations the approximate values p i substitute for 
exact values Pi. 
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Let mean number qk of particles of the fcth type (k = 1,2) 
that leave the pair of cells (1,1), (1,2) be called the particles of 
kth type flow intensity and mean number hk of particles of the 
kth type leaving the pair cells with a change of lane be called 
the changes intensity of the particles of the kth type. 

There is one particle, which can move, in the pair of cells 
(1,1) and (1,2), when the state of the four cells is one from the 
states £4, E w , £13, Eiq, En, E 2 g, E29, £30, £31 • If the state 
of the four cells is Eg then there are two particles which can 
move. Hence the value of q\ can be calculated as 

qi = (p4 + 2p 9 +Pi +Pl2+Pl3+Pl6+Pl7+P28+P29+P30+P3l)A i l- 

(6) 

Similarly we obtain the relations 

h = (pu +^13)^1, (7) 

Q2 = (P5+Pl0 + 2^i 1 +^i4+pi5+J?i 8 +^ig+p32+J?33+J?34+^35)/i2, 

(8) 

h 2 = (pu + P\b)P2- (9) 

Let q k , hk (k = 1, 2) be values of q^ and hk, which are obtained 
if in equations (6)-(9) pi substitutes for 

If iii > p,2 then q~i is some greater than the value q\ obtained 
by simulation and in contrast q\ is some smaller then the value 
(/2 obtained by simulation. As it be shown below the difference 
between the calculated approximate values and data obtained 
by simulation is smaller if values q\, hk determined by means of 
the correcting coefficient substitute for values qk, hk 

qk = (pk/pk)q~k, k = l,2, (10) 

h = (Pk/Pkjh, k = 1,2. (11) 

Let us explain why the correcting coefficient improves he ap- 
proximation. Let the ratio of the mean number of transitions of 
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particles of the kth type to number of the particles of this type 
be called the velocity of the particles of the kth type. Let the 
velocity of particles of the kth type be v k . Then 

v k = q k /{2p k ), A; = 1,2. (12) 

Because of the relation q k = 2p k -v k the obtained approximate 
value of the particles of the kth type differs from the exact value 
by the coefficient q\j p. Coefficient p/q 2 compensates the error. 

Suppose 

v k = q k /(2p k ), (13) 
v k = q k /(2p k ). (14) 

5. Approximations by means of Bernoulli scheme 

Let us consider another manner of evaluation of the steady-state 
probabilities and the macro-characteristics that is more simple 
than the manner described above. Assume that the probability 
of the state of each cell is independent of the states of other 
cells, and the cells with probability p\ contains a particle of the 
first type, with probability p 2 contains a particle of the second 
type and with probability 1 — p the cell is empty. Let pi (i = 
1,2,... 45), Vk, q k , h k (k = 1,2) be values of the appropriate 
characteristics that are obtained when these assumptions are 
used. 

Then p 1 = (1 - p) A , p 2 = p A = 2pi(l - p)'\ p 3 = p 5 = 
2p 2 (l - pf, Pq = h = pl(l - p) 2 , p 7 = pio = pi 3 = Pu = 

P17 = Pl8 = 2pip 2 (l - p) 2 , Ps = Pu = Pli 1 - P) 2 , P12 = Pl6 = 

2p 2 {l - p) 2 , p 15 = p 19 = 2p\{\ - p) 2 , p2Q = m = 2^?(i - p), 

P21 = P22 = P24 = P29 = P30 = P32 = 2/ofp 2 (l ~ p), P23 = P25 = 
P26 = P31 = P33 = PM = 2pip|(l - p), p 27 = P35 = 2^(1 - p), 
P36 = Pi, P37 = P39 = 2pfp2, ^38 = Pi3 = p\p\, PW = hi = ^p\p\, 
P\2 =PU = 1>P\P\, P45 = P\\ Vk = (1 + P~ p 2 )Pk, k = 1,2, 
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q k = 2p k (l - p)(l + p - p 2 )nk, k = 1, 2, hk = 2p k p(l - p) 2 pk, 
k = 1,2. 



6. Sistem solution and results of calculation 

Suppose pi = p2 = 0.25. 

The calculated approximate values are represented in table 
1. These values are compared to values obtained by simulation. 
The symbols for values which have been obtained by simulations 
are marked by star. The duration of the simulation interval is 
equal to 12000 time units, N = 500. 

We do not present the tables of calculations for p\ = 0.54, 
p 2 = 0.06 and p 1 = 0.08, p 2 = 0.72. 

The dependence of state probabilities on state numbers is 
represented on fig 2-10. 
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Figure 2: State probabilities for p\ = p 2 = 0.25, ^\ — 2, ji 2 — 1 

The calculated values of macro-characteristics are represented 
in table 2. 



7. Analysis of approximation accuracy 

Let us compare the values qk, vj-, hk, qk, Vk, hk, Qk: ^fc ? ^fc with 
the values obtained by simulation. 



Table 1: State probabilities, p\ = P2 = 0.25 
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Table 2: Macro-characteristics (the error that is relative to values obtained 
by simulation 
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q~k 


1.026(+19) 


0.058(-30) 


1.530(+15) 


0.057(-38) 


2.048(+26) 


0.056(-45) 


h* 


0.076 


0.003 


0.099 


0.002 


0.128 


0.001 


hk 


0.079(+4) 


0.004(+30) 


0.121(+22) 


0.003 (+48) 


0.159(+24) 


0.003(+220) 


hk 


0.207(+172) 


0.012(+284) 


0.311(+214) 


0.012 (+476) 


0.415(+224) 


0.012(+1052) 


hk 


0.082(+8) 


0.003 (0) 


0.126(+27) 


0.002 (0) 


0.168(+31) 


0.002(+100) 


* 

v k 


0.948 


0.692 


1.236 


0.767 


1.504 


0.842 


Vk 


0.915(-3) 


0.630(-8) 


1.356 (+10) 


0.704(-8) 


1.800(+20) 


0.747(-ll) 


Vk 


0.992(+4) 


0.496(-28) 


1.488(+20) 


0.496(-35) 


1.984(+32) 


0.496(-41) 


Vk 


0.950(+19) 


0.483(-30) 


1.417(+15) 


0.475(-38) 


1.896(+26) 


0.467(-45) 


pi = 0.08, p 2 = 0.72 


It 


0.044 


0.334 


0.039 


0.344 


0.048 


0.346 


qk 


0.045(+2) 


0.351(+5) 


0.050(+28) 


0.371(+8) 


0.056(+17) 


0.374(+8) 


qk 


0.074(+69) 


0.334(0) 


0.111(+186) 


0.334(-3) 


0.149(+210) 


0.334(-3) 


q~k 


0.072(+64) 


0.328(-2) 


0.105(+169) 


0.325(-5) 


0.140(+192) 


0.305(-12) 


h* 


0.002 


0.015 


0.003 


0.017 


0.004 


0.016 


hk 


0.004(+87) 


0.020(+36) 


0.004(+43) 


0.023(+34) 


0.006(+60) 


0.025(+54) 


hk 


0.010(+412) 


0.046(+207) 


0.015(+412) 


0.046(+171) 


0.021(+412) 


0.046(+188) 


hk 


0.006(+200) 


0.019(+27) 


0.009(+200) 


0.020(+18) 


0.016(+300) 


0.020(+25) 


«2 


0.275 


0.232 


0.244 


0.239 


0.300 


0.240 


Vk 


0.280 (+2) 


0.244 (+5) 


0.313(+28) 


0.258(+8) 


0.350(+17) 


0.260(+8) 


Vk 


0.464 (+69) 


0.232 (0) 


0.696(+186) 


0.232(-3) 


0.928(+210) 


0.232(-3) 


Vk 


0.450 (+64) 


0.228 (-2) 


0.656(+169) 


0.226(-5) 


0.875(+192) 


0.212(-12) 
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Figure 3: State probabilities for p\ = p 2 = 0.25, pi = 3, p 2 = 1 
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Figure 4: State probabilities for pi = p 2 = 0.25, pi = 4, p 2 = 1 



Suppose pi = P2 = 0.25. The values gi exceeds the values 
<7i (^i)- The relative error makes up the values from 19% to 70% 
on the three considered set of data. The value q~i (v\) exceeds 
the value q\ {v\). The difference q\ and q\ makes up from 37% 
to 116% of q\. The value q\ [v\) exceeds the value q\ (v*). The 
relative error make up from 45% to 113%. 

In the other cases the approximate values calculated accord- 
ing to formulas (10), (11), (14) are also nearer to value obtained 
by simulation than the values calculated by the two others man- 
ners, or in some cases the formulas (10), (11), (14) yields approx- 
imation only few worse. 

In each the three cases the value of q\ (£2) difference the values 
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Figure 5: State probabilities for pi = 0.54, p 2 = 0.06, ^\ — 2, // 2 = 1 
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Figure 6: State probabilities for p\ = 0.54, p 2 = 0.06, pi — 3, /i 2 = 1 

#2 (^2) * s n °t greater than 5% of q\. The value /ii difference from 
the value h\ not greater than 32% of h*. The value h\ exeeds 
h\. The relative error is not greater than 23%. 

Usually the relative error is the better the smaller is the ratio 
of jii to p 2 (/ii > ^2). 

The values pi, P2, ■ ■ ■ ,^45 that satisfies the system of linear 
equations usually differs from the values obtained simulation less 
than values calculated by means of Bernoulli scheme. 

The analysis of calculation for the cases p\ = 0.54, p<i = 0.06, 
and P2 = 0.06, p\ = 0.72, also shows that meaning Vk, (ft, hk, are 
usually more close differs from the values obtained by simulation 
less than values v k , q k , hk, v k , q k , h k . 



13 




1 2 3 4 5 6 7 8 9 1011 1213141516171819 2021 2223 24252627 2829 3031 3233 3435 3637 3839404142434445 



Figure 7: State probabilities for p x = 0.54, p 2 = 0.06, /ii = 4, /i 2 = 1 
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Figure 8: State probabilities for p\ = 0.08, p 2 = 0.72, \i\ — 2, // 2 = 1 



For /?i = 0.54, /92 = 0.06, the largest difference between q^ 
(vk) and gjj: (v^) corresponds to the case k = 1, fii = 4, fi2 = 1- 
In this case the value exceeds the value g|. The relative error 
is not greater than 20% . 

For Hi = 0.08, fi2 = 0.72 the largest difference between q~k 
(vk) ql (vl) is reached in the case k = 1, ji\ = 3, \±i = 1. In this 
case the value q\ {v\) exceeds values q{ {vl). The relative error 
is not greater than 28%. 
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Figure 9: State probabilities for p\ = 0.08, P2 = 0.72, p,\ — 3, //2 = 1 
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Figure 10: State probabilities for p\ = 0.08, p2 = 0.72, /ii = 4, /i 2 = 1 



8. Analysis of results of calculations of 
macro-characteristics 

If Bernoulli scheme approximation is used then calculated values 
of macro-characteristics are proportional to the particle transi- 
tions intensity. 

On the other hand it is clear intuitively that the fast parti- 
cles occur more often than slow ones in the states in that their 
transitions are prevented by other particles. Hence the ratio of 
velocity of particle to its transition intensity (fit) for the fast 
particles is less than for the slow particles. It is confirmed both 
by simulation and calculations. The difference between v<ij \xi 
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and Vi/ni is the greater the greater the /ii//i2 is. Because of 
that for fixed fi2 both value v* and value v\ calculated accord- 
ing to formulas (10), (14) increases slower than linearly if /ii 
increases. But it can be noticed that for the valued obtained by 
simulation the slowing down of increasing of the velocity of the 
particle of the first type with increasing of fi\ is more obvious. 
It can be noticed also if \i<i is fixed and fi\ increases then the 
values V2 and v\ increases slowly and the value V2 is constant. 
For fixed \i<i and increasing fii the values h\ and h\ increase ap- 
proximately linearly. The ratio of hk/q~k to h%/q% for the fast 
particles is greater than for slow ones (moving forward the fast 
particles change lane more often than slow particles). 

9. Conclusion 

It is presented a model of the traffic flow on two lanes. An ap- 
proximate manner is elaborated for calculation the steady state 
probabilities and macro-characteristics of model. The results of 
calculations by means of the elaborated method are compared 
with results of simulation. This comparison shows that the pre- 
sented method yield rather good approximation and describes 
the qualitative behavior of the considered characteristics. 
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